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1 Introduction 

The non-perturbative aspects of the strongly interacting supersymmetric 
(SUSY) gauge theories were intensely investigated in the past [|l|] and re- 
cently they have been the object of renewed interest @, 0]. These theories are 
interesting from a phenomenological point of view as their non-perturbative 
properties might play a crucial role in the understanding of the SUSY break- 
ing mechanism |l| . However, besides the relevant phenomenological implica- 
tions, the SUSY gauge theories have intrinsic importance as their very nature 
allows the calculation of some exact non-perturbative solutions [^ ^ . 

In this paper we concentrate our attention on the most simple SUSY 
gauge theory: the pure A^ = 1 SUSY Yang-Mills (SYM) with SU{N^) gauge 
group. This contains the purely gluonic action, plus one flavour of Majorana 
fermions in the adjoint representation of the colour group. It is believed 
that this theory is characterised by the same non-perturbative phenomena 
as QCD: colour confinement and chiral symmetry breaking |^. Nevertheless, 
noticeable differences between the SYM theory and QCD appear to exist 
even at the fundamental level. Due to supersymmetry there is indeed a new 
anomalous SUSY current which belongs to the same supermultiplet together 
with the anomalous chiral and energy- momentum tensor currents 0|. 

The fundamental question of the breaking of the supersymmetry in A^ = 1 
SYM theory was addressed in 0, 0|. According to the general argument 
of the Witten index f^ or the Veneziano-Yankielowicz low energy effective 
theory |Q, the supersymmetry does not break. Nevertheless, here the chiral 
symmetry breaks and the gluino condensate acquires a vacuum expectation 
value [^, 0. However, the authors of Ref. [0 argue that in this theory the 
spontaneous breaking of chiral symmetry implies the spontaneous breaking 



of supersymmetry due to non-perturbative effects. 

From 0, ^ we learn that the low energy supermultiplet contains three 
degenerate massive colourless bound states: a scalar, a pseudoscalar and a 
fermion field (where the appearance of the fermion field in the low-energy 
supermultiplet is a consequence of the colour adjoint representation for the 
gluino). Moreover, there is no pseudo-goldstone boson (or pion) associated 
with the chiral symmetry breaking, as the latter is broken by the anomaly. 

Some time ago Curci and Veneziano [§] suggested that the SYM theories 
can be studied non-perturbatively on the lattice by using numerical Monte 
Carlo simulations. This is clearly analogous to the approach used in QCD 
theory. They argued that, even if the lattice breaks explicitly supersym- 
metry, it is possible to recover the SUSY and chiral Ward identities on the 
continuum limit. 

Recently two different collaborations 0-|T^ studied non-perturbatively 
on the lattice the spectrum of the SYM (A^ = 1) theory following the guide- 



lines suggested in ||g]. In [|Tl|], because of the limitations deriving from the 
use of computing resources, the quenched approximation was used to study 
the spectrum. This approximation consists in neglecting the internal gluino 
loops. In other words, in the correlation functions of fermion composite op- 
erators, the fermion determinant is not included. In SYM theories, if general 
arguments are taken into account [^], the quenched approximation cannot 
be justified on the basis of large A''^ dominance, since gluinos are in the ad- 



joint representation of the colour group. However, in [|T^ it is observed that, 
within the statistical errors of the Monte Carlo simulation, the spectrum 
seems to show no deviations from the supersymmetry expectations under 



the approximations considered there. In connection with this result, in |jT3[, 



by means of an effective lagrangian approach, the splitting in the low energy 



supermultiplet induced by the quenched approximation has been analysed, 
and it is argued that the sphttings from the supersymmetry predictions are 
smaU. 

In any case, in [jlO|, |l^ unquenched results for an SU(2) SYM gauge theory 



were obtained by using numerical simulations with dynamical gluinos. It is 
likely that new results in this field will quickly follow. 

In this paper we study the spectrum of the SYM (A^ = 1) on the lattice 
at strong coupling and in the large Nc limit. The lattice strong coupling ex- 



pansion technique (see [Ti\ and references therein) has been extensively used 
as an analytical probe to test quahtative properties both of the continuum 
and of the lattice theory by itself. Up to now, however, the main part of 
the results refers to QCD, and -to our knowledge- no work deals with a su- 
persymmetrised version of it. In detailed computations the strong coupling 
expansion of SU{Nc) theories is frequently combined with the large Nc one. 
Our reference comes from works in which issues such as the phase struc- 
ture of QCD or the computation of meson and baryon masses are addressed 



15[-|21|. The most frequent computational frameworks split into two types: 



Effective actions: The Wilson-Dirac lattice action is considered at large 
Nc and small f3. The large Nc expansion can be recognised as a saddle- 
point expansion of the gauge functional integral, previously simplified 



by the (3^0 limit |T5|-[0. The method allows the study of the large 
and small hopping parameter regions at the expense of introducing 
assumptions on the form of the saddle points. The main disadvantage 
of this method for our purposes is that the one matrix integration 
formulae, a basic ingredient of the construction, are not available for 
matrices in the adjoint representation, and the generalization does not 
seem straightforward. 



Path resummation: The ferinion matrix M is inverted by using the 
standard hopping parameter expansion M, |T6|, |21|. This allows to 



compute propagators M"-*- in terms of sums over paths on the lattice 
(random walks); the objects to be summed are traces of products of 
spin matrices (rl ± 7^), r being the Wilson parameter. We will see 
that provided the hopping parameter k is small enough and some con- 
straints on the parameters are imposed, no assumptions are required 
to perform this computation in the N^ -^ 00, (3 -^ Q limit. In some 
of the previous references the analysis has been performed only for the 
case r = 1, where the structure of the series is considerably simplified 



16]. The main difficulty of path combinatorics arising for r 7^ 1 was 



addressed in Ref. |^. In Ref. [^ an independent derivation of the 



main resummation formulas for r 7^ 1 is given. 

In our case, it will be shown that the hopping parameter expansion in 
terms of random walks is valid, with slight changes, for the case of gauge 



fields in the adjoint representation. Thus, using the formulas of Ref. |23| we 
will be able to give expressions for propagators and masses for any value of 
r and small enough k. We considered that keeping r arbitrary could be very 
important. This allows the possibility of searching for multicritical points, 
where several lattice masses vanish. In particular, we investigated whether 
one could find a critical point in the k — r plane where the massless modes 
would form a supermultiplet. This would signal a possible candidate for a 
supersymmetric continuum limit, in the spirit of the chiral restoration of or- 
dinary QCD. For example, if a low energy theory of the type described in ^ 
would take place, one would find a massless pseudoscalar meson, together 
with a scalar meson and a fermion. We should warn the reader that by 
massless we refer to the vanishing of the lattice masses, not necessarily the 



renorinalised physical masses. However, in our approximations the diagrams 
contributing to the anomaly, and giving mass to the pseudoscalar, are sub- 
leading, indicating that the physical masses are also zero, as predicted for a 
Goldstone boson associated to the spontaneous breaking of chiral symmetry. 
This situation contrasts with what one obtains in perturbation theory and 
constitutes one of the most salient features of our result. 

In addition to the path resummation formulae, we also need to analyse 
the behaviour in the large Nc limit of integrals over the group of products 
of matrices in the adjoint representation. These are studied in |^. In this 
paper we will only need to know the order in N^ of certain integrals. Given 
this, the actual leading order results are very simple. 

To increase the usefulness of our paper, many of our formulas will be 
given for arbitrary space-time dimensionality d. In addition, we will indicate 
the necessary changes to make the formulae valid for matter fields in the 
fundamental representation and for Dirac, rather than Majorana, fermions. 
Specific attention will be paid, however, to the supersymmetric cases oi d = 
3,4. 

The paper is organised as follows: in section 2 we present the formalism 
and general formulae for the expectation values and correlation functions 
of two and three gluino operators in the strong coupling large Nc limit. In 
section 3 we analyse the specific channels and present the results for the 
spectrum in the meson and fermion sector. Finally in section 4 we summarise 
our conclusions, discuss the physical implications and explain the prospects 
for possible future extensions of our results. The paper is completed by 
collecting in the Appendix the terminology and formulae on lattice paths 
that we will use along the text. 



2 General formalism 

We begin by fixing our notation and terminology. Let us consider a d- 
dimensional hypercubic lattice C = Z'^. We introduce an index set / having 
2d elements. To any direction in space-time we associate two indices /i and Jx 
corresponding to the two opposite senses (forward and backward) associated 
to each direction /i. Given an element a G /, the element a labels the one 
with reverse orientation {fl = fi). To any element a G /, we can associate a 
vector V{a) as follows: 

V(/i) = e^ Vm) = -e, , (1) 

with e^ the unit vector in the fi direction (the lattice spacing is set to 1). In 
the Appendix we give some additional results and terminology that we will 
be using in the following sections. 

Let us now write down the lattice version of the SUSY Yang-Mills action: 

S = pS,+ ^-^,^,M,, , (2) 

where (3Sg is the pure gauge part and \E'j is a Grassman variable representing 
the field of a Majorana fermion. The index i is a short form for the three 
indices n, a and A. The index n specifies a point in the d-dimensional 
hypercubic lattice C. The index a takes A*"^ — 1 values corresponding to the 
dimension of the Lie Algebra of SU{Nc). The index A is a spinorial index 
taking 2['^/^l values (the symbol [x] stands for integer part of x). Without 
much problem, but at the expense of breaking supersymmetry, we could add 
a flavour index ranging over a finite number of values. The matrix M must 
be antisymmetric and its form is given by 

M = C M (3) 

M={l-Eaei^^) , (4) 

7 



where I is the unit matrix and Cab is the charge conjugation matrix, satis- 
fying: 

7*C = -C7, (5) 

C* = -C , (6) 

where the superscript t stands for transpose. 
Finally the matrix A^ is given by: 

(Aa)ij = K5mn+V(a) U^i^) ('"I " la)AB , (7) 

with i = {n,a,A), j = {in,h,B)\ k is the hopping parameter, and r is the 
Wilson parameter. For a = /i, U^{n) is the gauge field link variable (which 
is in the adjoint representation) and 7^ the Dirac matrix, while for a = /i we 
have Ufi{n) = Uj^^n — e^) and 7^ = — 7/^- Notice that we have: 

A^A^ = K\r^-l)I . (8) 

It is easy to show that given the conditions Eqs. (^)-(^), the matrix M is 
indeed antisymmetric. 

Customarily, the value of r is taken in the interval [0, 1]. This follows from 
the requirement of Osterwalder-Schrader positivity. However, at strong cou- 
pling this is not a necessary condition for the physical correlation functions 
to admit an analytical continuation to unitary Minkowski Green functions. 
For that reason we will be working for arbitrary r and comment at the rel- 
evant places whether the p-gluino correlation functions satisfy the positivity 
conditions. 

The constraint imposed by the existence of the matrix C, together with 
the needed matching of fermionic and bosonic degrees of freedom, makes 
the naive continuum limit of the above lagrangian supersymmetric (A^ = 

8 



1) in d = 3,4. The same lagrangian is supersyinmetric in d = 10 if \l/ 
is a Majorana-Weyl field, and in (i = 6 if \1/ is a Weyl spinor. However, 
the requirement of Weyl character conflicts with the well-known absence of 
chirality on the lattice, which prevents us from directly writing a lattice 
version of the required lagrangian (indeed, our Wilson-type action breaks 
chirality explicitly). On the other hand, the general argument by Curci and 
Veneziano linking the recovery of chirality and supersymmetry could play 
a role in the interpretation of these cases, as of course does in the (more 
transparent) ci = 3, 4 as well. 

Having in mind this caveats regarding the interpretation of the results, in 
the rest of our exposition we will try to work as much as possible in arbitrary 
dimension d without specifying more. This has the advantage that it allows 
our formulae to be used with little changes by other researchers with other 
physical interests. For example, links can be established with the technique 
of 1/d expansions |^ . 

We will concentrate upon gauge invariant operators of the form: 

0,(a;) = v]/- (x) . . . ^Xix) Sf^-^^ C^^ , (9) 

where C"'^"'"'^ is an invariant colour tensor and S^ ^"' ^ a spin tensor. The 
index i labels different possibilities for the definition. The restriction to 
operators obtained by multiplying gluino fields at the same lattice point (ul- 
tralocal operators) will not affect our spectrum results which will be general. 
In particular, this affects the composite fermion operator which belongs to 
the same supermultiplet as the scalar and pseudoscalar: 

x^(^) = ^K.)^^^;.(^)n(^) , (10) 

where in this formula the fermion fields x 5 ^ ^ind the Yang-Mills field F^^/ 
live in the continuum spacetime, and cr^u = \['y^,'yu]- Later on, we will 

9 



argue that its corresponding minimal mass will be contained within the set 
of masses corresponding to 3-gluino operators. 

Note that the operators in Eq. |^ are non- vanishing only when the gluino 
fields are combined in a completely antisymmetric way. By putting this 
together with the requirement of Cj being an invariant tensor, we arrive at 
certain constraints on the possible operators; for example, in the case of 
meson-like operators {p = 2) it turns out that, necessarily, Cj = I and Si is 
an antisymmetric matrix. As an example. Table 1 displays a complete basis 
of 2-gluino operators in d = 4. 

In this paper we will be interested in computing the expectation values 
of these operators and products of these operators: 

{0.{x)) = ^l[( fdnn)l[fdU,{n)\o,{x)e-' (11) 

neC \-^ fJ. -^ J 

G,,{x-y) ^ (0.(x)0,(i/)) = 

i n (/^^H UJ^^^^'^yj 0.{x)0,{y)e-' (12) 

ne£ \ M / 

at strong coupling. We will be able to accomplish this goal for /5 = 0, and in 
the large Nc limit. Corrections to the formulae in powers of /? and -^ are in 
principle feasible and will be considered elsewhere. 

As usual, one can explicitly integrate out the fermions in Eq. (p!TD -([l2D. 
The main formula that one uses is: 

(13) 
where P/(M) stands for the Pfaffian of the matrix M. The square of the 
Pfaffian is the determinant, so that up to a sign we can write: 

P/(M) = y/det{C)det{M) = exp{^Tr(log(A^))} , (14) 

10 



where we have used the standard exponential representation of a determinant, 
and the fact that we can choose ^/det{C) = 1. Concerning the sign, it is 
clear that what matters is the relative sign for different values of the gauge 
field configuration. However, for very small k, the matrix A4 is close to the 
unit matrix and our representation (Eq. (P^)) is valid. Problems can only 
arise when one of the eigenvalues of Ai crosses zero. Using the Schwartz 
inequality and our expression for A4, one can conclude that this problem 
never occurs provided \k\ < 2d(\r\+i) - 

What we will do now is to expand the quantities entering in Eq. (|1^) as 
a sum of paths. Using the terminology of the Appendix, we obtain: 

{M-\x, y))% = E,esi.-,y) W'^'i^) ^ABil) (15) 

P/(M) = exp{i E.e£ Er=i E,.5,(.^.) i Tr{Wm Tr{T{j))} , (16) 

where x,y are lattice points, W{'y) is the path ordered product (along the 
path 7) of the gauge field link variables, and r(7) denotes the appropriate 
product of the spin matrices: 

r(7 = {x, a)) = K\r - 7,J . . . (r - 7„J . (17) 

Let us first consider 2-body operators: 

Oi(x) = Sf^ ^^(x) ^%{x) = ^(x)^i^(x) (18) 

with S^ = C-^Si . (19) 



The second expression in ([T8|) has the same form as for Dirac fermions, the 
difference being that for Majoranas we have the relation \E' = \E'*C. Notice 
that Si is an antisymmetric matrix. Using the formula (O) and expanding 
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by using Wick's theorem we obtain: 

(0,W) = - ' {Tr{M-\x,x)C-^S,) P/(M)>, (20) 

[-2 Tr [C-^ S,M'\x, y) C"^ SjM-^{y, x)) (21) 

+Tr {M-\x,x)C-^ Si)Tr {M~\y,y)C-^ S,)])g , 

where Tr denotes the trace over colour and spin indices and the symbol 
( )g denotes expectation value with respect to the pure gauge action part. 
The first term on the right hand side of Eq. (^) represents the so-called 
OZI contribution, while the second term contains the disconnected and non- 
OZI contributions 0. Now using the formulae (|15D-(|16D we can write back 
expressions ([20|)-(^l|) in terms of products of sums of paths. 

Up to now the expressions have been completely general, but now we 
will consider the simplification arising from considering /5 = and N^ ^ oo. 
What we need to know is the expectation value of the product of traces of 
Wilson loops in this limit. The main result that we will use |2^ is that for 



a collection of closed paths {71, . . . ,7s} which are not pure spikes (see the 
Appendix) we have: 

(rr(H^(7i))...Tr(H^(%))). = 0(l). (22) 

This is true both for connected and disconnected expectation values. How- 
ever, for a pure spike path 1^(7) = I, and therefore its trace is A^^ — 1. We 
see that the connected correlation between a pure spike Wilson loop and any 
other operator vanishes. Our first conclusion is then that, at leading order 
in Nc, the factor P/(M) cancels between numerator and denominator. This 
is precisely the quenched approximation, which turns out to be exact in this 

12 



limit. In principle, the result is surprising since fermions carry indices taking 
A*"^ — 1 values (as gluons), and the usual arguments why fermion loops are 
subleading in the fundamental representation do not apply here. This adds 
to the results obtained by other methods [|ll|, |l^, pointing towards the fact 
that the deviations introduced by the quenched approximation are not too 
large. 

Other conclusion that follows from formula ( ]22|) is the suppression in the 
large Nc limit of the non-OZI contributions to the connected correlation func- 
tions. This includes the mass induced by the anomaly on the pseudoscalar 
Goldstone boson. We should, hence, expect such a massless state signaling 
the recovery of chiral symmetry and its spontaneous breaking. As mentioned 
in the Introduction this result, as well as the exactness of the quenched ap- 
proximation, contrasts with the behaviour obtained in perturbation theory 
to leading order in Nc- 

Our results have been obtained at /? = 0. However, one can introduce the 
pure gauge action in various ways having the same naive continuum limit. If 
one chooses to write the Wilson action in the adjoint representation of the 
group, formula (^) implies that to the order we are working (and in the 
small /3 region) there are no corrections to any order in (3. This is not the 
case if the customary fundamental representation Wilson action is chosen. 

The conclusion of the previous paragraphs is that to leading order in A''^ 
all that we have to take into account are closed loops which are pure spikes. 
For that purpose the results of Ref.||23[ and collected in the Appendix are 
needed. Notice that for a pure spike path of length L, the spinor matrix is 
given by: 

r(7) = (fi:2(r2-l))^/2l. (23) 
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Thus, using the formulae given in the Appendix we can conclude: 

(Oi(x)) = -{Nl-l)F{Q,K\r^-l))Tr{C-'Si) = 



1 



-[N',-l)Tr{S.)- , (24) 



_ 2d c 

with 



This is the contribution of closed paths which are pure spikes. The corrections 
coming from other paths are order 1, and thus subleading. The previous 
formulae are obtained by resummation of a series. The radius of convergence 
is given by the closest singularity. Thus the formulae are strictly speaking 
only valid in the region [/{^(r^ — 1)| < ^,^j_-^, . It is possible that, by analytical 
continuation, the formulae could be valid in some points beyond this region, 
such as for larger negative ^^(r^ — 1). 

Expression (|2^ ) is also valid for Dirac fermions. It is also valid if the 
fields (Dirac or Majorana) transform in the fundamental representation of 
the colour group, provided {N^ — 1) is replaced by Nc, the dimension of the 
representation in question. 

Now we look at the correlation function of two-fermion operators. In this 
case we have two factors of M. and hence we have an expansion in terms of 
paths 7 going from x to y, and paths 7' going from y to x. Nevertheless, 
the integration over the gauge group forces the combined path to be a pure 
spike path. To take this into account we proceed as follows. We replace the 
summation over paths by a summation over paths with no spikes, resumming 
all paths which have such a path as reduced path. Thus, each term in the 
new expansion corresponds to a reduced path 7 going from x to y, and 
another one 7' that returns to x. However, now the condition imposed by 
the integration over the gauge group is simply that 7' is the reverse path of 7 

14 



(which we label 7~^). In this way the double summation reduces to a single 
summation. Summing up all that we have just expressed in words, we can 
give the following formula for the connected correlation function: 

oo 

Tr(s,T{^)S,T{r')) . (26) 

The previous formula has been written in a way which makes it valid for 
Majorana [r] = 2) or Dirac [r] = 1) fermions. The symbol Dpt stands for the 
dimension of the gauge group representation (A^^ — 1 for the adjoint and Nc 
for the fundamental). Now we can use the expression for F{L,K,'^{r'^ — 1)) 
and the formulae for resumming over paths that are given in the Appendix, 
to conclude: 

(27) 



where we have: 



1-M4a: 



R^ix) ^ -vDn- ^^ (28) 

1 — T 

e.W ^ (i-r^+p^ (29) 

M'f) = K'5Ze*"(r-7„)®(r-7„) (30) 

C2 = C®C. (31) 

\Si) is just given by the matrix iSj, but considered as a 2'^^^ ■ 2^2^ dimensional 
vector. It is useful to express the matrix elements of the 2^[2l x 2^t2l matrix 
A2 between Si states in terms of the matrices Si defined in Eq. (^). We 
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have: 

(5,1 k^iif) C2 1 \S,) = i^'Yl ^""^ ^^['^^(^ - ^°)'5'^-(^ + ^-)] • (32) 
With this interpretation, formula ( pT]) is vahd for Dirac fermions as welL We 



will leave to the next section the evaluation of this expression and the study 
of the properties of the resulting propagator. 

We now turn our attention to 3-gluino fermion operators of the form: 

In this case there are two possible invariant colour tensors: dabc and fabc- 
The main result that we will need on the group integration at large N^. is the 
following: given three paths without spikes 71, 72 and 73, we have: 

+subleading terms , (34) 

where C"'^'^ is either f or d { the antisymmetric and symmetric SU{Nc) invari- 
ant tensors). The mixed terms f — d are sub leading in A^^- Also subleading 
are contributions in which the three paths are non-equal. Using this expres- 
sion and the formulae derived in the Appendix, it is possible to compute the 
expectation values of products of 3-gluino operators. Now the Wick expan- 
sion gives a total of 6 terms (once the A'"c-subleading ones are discarded). 
These terms ensure that if the colour invariant tensor is / or d, the spin ma- 
trix Sj can be chosen completely symmetric or antisymmetric respectively, 
as required by Fermi statistics. With this choice the 6 terms give the same 
contribution, and the correlation function can be written as: 

(0.(x)0,(y)) = RsiO n( / ^) ^^"^^"'^ ('^.l [93(01 - A3(^)]-^C3 ^ \S,) , 

16 



(35) 



where now: 



, (l-x)3-(-£^)3 

^3(.) . -^^' \,_\^^~y! (36) 

Qsix) ^ (l-^)'+(^^3Ij^ (37) 

A3((^) = «:'5^e*^'^(r-7,)®(r-7j®(r-7,) (38) 

Cs = C(^C^C. (39) 

As before, the vector \Si) is the one constructed from the corresponding spin 
matrix. Exphcitly: 

{S,\AMC,'\S,) = ^^Eaei^'''"^ (40) 

Sf'^'^^'iir - 7„)C-i)^i^i((r - ia)C'Y''''{{r - 7a)C-i)^^'^^'5f ^^^^ . 

The class of operators considered (Eq. (|33D) does not include the lattice 
counterpart of that in Eq. ([T0|). A possible candidate would be: 



X^(x) = l{cr,.)^''v;t{x)r'^mx) , (41) 

where Vf^^ix) is an appropriate combination of adjoint plaquettes in the 
(yU, z/) plane whose naive continuum limit is, up to a convenient multiplicative 
factor, the adjoint gauge field F^i,{x). This has the advantage of including 
only gauge variables in the adjoint representation, allowing the use of our 
integration formulas. However, by examining the strong coupling large Nc 
expansion of the correlation of 2 such operators, one easily realises that it is 
given by combinations of triple paths joining the two operators. Thus, we 
expect that the mass spectrum following from 3-gluino ultralocal operators 
would include also the states coupled to (|4lD . 
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To conclude we simply want to mention that in a similar way one can 
obtain expressions for expectation values and correlations of p-gluino opera- 
tors. No additional difficulty arises, and the final expression looks just like 



Eq. (pSD but with corresponding functions and matrices -Rp, Gp, Cp, and Ap. 



In particular: 



e.W - ('--r + j^j^j^ (42) 

A,{^) = k' Ve*-(r-iJ»...®(r-7j . (43) 
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3 Explicit results on the propagators and spec- 
tra 

In this section we will analyse the results on the expectation values presented 
in the previous one. Our main goal will be the extraction of the spectrum of 
the theory at /? = and leading order in the 1/Nc expansion. 

We will first of all look at the expectation values of single 2-gluino oper- 
ators. Our main result is formula (0). The only independent operator Si 
giving a non- vanishing spinorial trace can be chosen as Si = CI. On physi- 
cal terms, it corresponds to a non- vanishing gluino scalar condensate for the 
full range of values where our resummation is valid (^—3 — ^ — I)' O^^e 
the different normalisations of the fields and operators, and the appropriate 
colour factors are taken into account, this expression coincides with the re- 



sult given in Ref. ^^ for this expectation value when the fermion field is in 
the fundamental representation. 

Our next step will be to analyse the results for the correlations of two 
fermion operators Gij{x), and the corresponding {meson) spectrum. The 
expression given in the preceding section for the correlation (formula (p7|)) 
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requires the inversion of a 2*^ x 2'^ matrix. (For odd space-time dimensions one 
must replace d by 2[|]). This matrix is Q2{C)^-^2{'^) defined in Eqs. (p9|,pOD. 
Thus, in even space-time dimensions, it is convenient for the study to choose 
as a basis of the 2°'-dimensional space of meson operators, those correspond- 
ing to Si being the standard basis of the d-dimensional Chfford Algebra, for 
which we adopt the writing: 

^K) = e^'("'')7o°---7."-I^ , (44) 

where 5{n^) is an appropriate phase which we will choose equal to zero in 
what follows. Thus, one state of the basis is the scalar operator Ss corre- 
sponding to the unit matrix, other elements are the vector operators Sy(^) 
corresponding to the gamma matrices 7^, and so on. 

For odd space-time dimensions, one could also consider the operators 
S{n^) associated to the standard basis of the Clifford algebra (which is basi- 
cally the same as for dimension d — 1), but one must take into account that 
they are not independent. They are constrained by the identity: 

7o---7d-i = KdI , (45) 

where K^ is a phase depending on the space-time dimension. 

Going back to the even- dimensional case, we have to express the matrix 
elements of the matrix A2(v?) within this basis. To do so it is convenient to 
view the 2'^-dimensional space in question as the Fock space of a system of 
fermions: the gamma- fermions. Each integer n^ G {0, 1} entering Eq. (^^ 
can be interpreted as the occupation number of the state /x. It is convenient 
to add an additional one-particle state labelled '—1' whose usefulness will be 
clear in what follows. Thus, in the standard second quantization notation 
we can write: 

S'(n^) = |n_i,no,...nd_i) , (46) 
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with S{n^) defined in (0). The extra occupation number n_i is fixed to be a 
parity bit state, taking the value 1 when the total number of gamma-fermions 
in the other states is odd and the value if it is even. Hence, in both cases, 
this imposes the constraint that the total number of fermions must be even. 
With this convention it is possible to express the matrix A2{(p) in terms of 
creation and annihilation operators of these fermions as follows: 



A^i^) = K\2d{r^ - 1 + 2a+ia_i) - 2(4a+ia_i - 2) E;.=oCos(^^)a+a^ - 

-4^'"E^=isin(^^)(a+a_i + a+ia^)) , (47) 

with a = Yfi7=o cos(<^^) . 



We see that the operator conserves the number of gamma-fermions. Hence, 
each even number of gamma-fermions 2p characterises a block in which A2(v5) 
can be diagonalised or inverted. Within each block one has two subspaces 
corresponding to n„i = 1 and n_i = 0, which correspond to the product of 
2p — 1 and 2p gamma matrices respectively. The matrix A2{(f) mixes these 
2 subspaces. 

Let us clarify the previous formulae by looking at a few examples. Up = 0, 
we have the gamma-fermionic Fock vacuum state, which is an eigenstate of 
A2(v5) with eigenvalue 2aK^(r^ ~^)- This state is precisely the scalar meson 
operator. Next, we consider the space of 2p = 2 gamma-fermions. The 
subsector n_.i = 1 corresponds to the vector operators (5'y(^) = 7^) and the 
n_i = to tensor states {St{pu) = 1/yu)- The conclusion is that vector and 
tensor states mix between themselves but not with other states. Considering 
the space of 2p gamma-fermions, we conclude that the operators associated 
to the product of 2p — 1 Dirac gamma matrices mix with those involving 2p 
gamma matrices, but with no other states. The inversion or diagonalisation 
problem has been considerably simplified with this technique. This is the 
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generalisation of the block structure found by previous authors studying QCD 
at strong coupling |16[ . 



In the previous analysis we have not taken into account the restriction 
imposed by the fact that our gluinos are Majorana. As mentioned before, in 
this case if the operators involve the product of the gluino fields at the same 
point (ultralocal operators), the matrices Si can be chosen antisymmetric. 
This restriction translates in our language into p being an even number: the 
number of gamma fermions must be a multiple of 4. Thus, for instance, 
the only relevant blocks in c? = 4 for this ultralocal operators are the scalar 
singlet and the p = 2 containing the pseudoscalar (75) and the axial vector 

We now introduce an important symmetry of the operator given in (^7|). 
This is the unitary transformation C related to the charge conjugation of 
gamma- fermions: 

Ca+Ct = a^ (48) 

Ca+iCt = -a_i (49) 

C|0) = K^|1,...,1) . (50) 

Up to a phase the operation exchanges occupied by empty for all states. One 
can easily see that the previous transformation commutes with the operator 



A2{ip) given in expression (^7]) . Notice, that with this change a state with 2p 
gamma- fermions changes into one with {d+l—2p), which for even space-time 
dimensions {d = 2s) is an odd number. We see that the spaces with an odd 
number of gamma fermions are useful after all. Thus, there is a hierarchy 
of complexity in the Fock space of gamma-fermions as the number grows. 
The Fock vacuum corresponding to the scalar operator is an eigenstate of 
A2{if). Next, comes the 1-particle space, which through charge conjugation 
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of gannna-fermions corresponds to the operator involving the matrix 7 = 
7o ■ ■ -7^-1 and 77^. In what follows we will proceed to invert the matrix 
02(01 ~ A2(v5) and obtain the propagator for these simplest cases. 

Before proceeding to the inversion, let us comment about the necessary 
changes to deal with an odd-space time dimension d. In this case we might 



introduce gamma-fermions as well, but due to the constraint (^Sf ) there are 
actually 2 states corresponding to the same operator. However, with a bit 
of effort one can show that the two states are precisely the 2 states that 
are mapped by the transformation C, provided K'^ is chosen equal to K^ 
entering in Eq. (53). With this choice, one can see that for odd space-time 



dimension expression (|47|) remains valid, but that the physical space of meson 
operators has to be identified with the subspace of the gamma-fermion Fock 
space which is invariant under C and has an even number of gamma-fermions. 
With this in mind all that follows can be applied to even and odd space-time 
dimensions equally. 

Now we consider the scalar state (corresponding to the unit matrix) first. 
Indeed, to comply with the normalisation chosen for the d = A case in Ta- 
ble 1, we take Ss = 2~2[2]l. The matrix A2(v?) reduces here to the constant 
2cr/t^(r^ — 1). Then we can directly write the momentum-space propagator 
explicitly: 



Gss{^) = ^ ... ^^1 — r^ ' (51) 



<^2(o-E';icos(^ 

where we have defined the following functions: 



-7] DR(2d-l) (1 + l^x) 
2x(l-x)(l-^x) 

^ ' 2x(l-x) 

and rj, ^ and Dpt are the ones defined in section 2. The function $2 is 
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decreasing for all ^ 7^ in the convergence interval, and in addition satisfies 
the following properties: 

$2(1/2) = d 

e e (0,1/2)^ $2(0 >rf 

^ < -^ $2(0 < -d . 

Next, for even space-time dimension, we proceed to study the space of 1 
gamma-fermion. As mentioned previously it corresponds to the matrix 7 (75 
in 4 dimensions) and 77^. Expression (^Tf ) reduces in this case to: 

d-l 

A2{<f) = 2K\r^ - l)a + K^ J2 T^pa+ap , (53) 



a,/3= 



where 



T_i_i=4a 
T_i^ = T^_i = Airsm{ip^) (54) 

Now the propagator can be obtained by inverting the matrix 62! — T. This 
can be done by making a non-unitary change of variables which brings T to a 
diagonal matrix up to a 2 x 2 block. In this way, one obtains the expression for 
the propagator in momentum space G^^"^-* in the 1-gamma fermion sector, i.e., 
the axial vector and the pseudoscalar block in d = 4. To make a contact with 
the usual conventions, we change the normalisation of the operators to the 
forms (remember d is even) P = 2~^^^'j and A(p) = i2^^^^'jp'j, generalising 
again the 4-dimensional ones given in Table 1. The result is: 



Gpp{^) = H{0{a, + J2- 



2^ 1 
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GpA{p){^) = Ga{p)p{^) = -—Gpp{ip) 



Otp 



Ga{p)a{p){v) = —[ad+ Yl —]Gppi(fi) 



"^ V ^=0,^^p "'^ 



GA{p)A(a){v) = GA(a)Aip){v) = ^-^Gpp{lp) , pj^a , (55) 

where H and $2 are the quantities defined in Eqs. (|52|), and the functions 
ad, dp, and /?^ have the following expression: 

r^ + 1 '^'^ 
ad = $2(0-^:^3^X1^°^^'^'^) (^^) 

p=0 

d-i 2 

ap = $2(0 - Yl cos(^p) - ^;^— Y cos((/?^) 

p=0 

Pp = 2^ — rsin((^^) . 
r^ — 1 

Obtaining momentum-space propagators for other blocks is feasible, but the 
expressions become more and more complicated. Furthermore, for ultralo- 
cal operators in 4 dimensions, the previous propagators are the only non- 
vanishing ones. Thus, we will focus in what follows on the analysis of the 
meson spectrum. 

The lattice masses are the minima of the lattice energies as we vary the 
spatial momentum 0. These minima can only occur at special momenta 
= (^{special) ^^{speaa } _ g^^^_ rpj_^g advantage is that now sm{(pf^^'^^''') = 



0, which simplifies expression (|47|) considerably. The procedure to obtain 
the masses is the following: Extract the eigenvalues of the matrix 02(^1 ~ 
A2(<^^*''*^'^*"'^), which are functions of the temporal momentum (^q. Then 
determine (^0° ^5 the (complex) value of (^0 for which the eigenvalue vanishes. 
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The lattice masses are now given by M = —log{\e'^'^o |). This coincides with 
the definition of mass as the exponent controlhng the decay of correlation 
functions at long times. 

Now let us proceed to obtain the eigenvalues. By looking at expres- 
sion (|47|) , one sees that the matrix is diagonal except for the term propor- 
tional to {aQ a-i + a^itto) . Thus, the occupation numbers of the spatial states 
Hi are not changed by the operator. Hence, for fixed values of these num- 
bers, the operator reduces to a 4 by 4 matrix: an operator acting on the 
two-state fermion system labelled by n^i and uq. Furthermore, the states 
having n^i = uq = and ra_i = rio = 1 are eigenstates of the matrix 
A2((^*-^^'^'^*"^^). The other two states are mixed, but finding the eigenvalues 
and eigenvectors is trivial, since it is a 2 x 2 matrix. We can summarise the 
results obtained in the following formulae: 

^^(^-(specmO)|o^O,n) = (2K2(r2-l) cos((/?o) + 2kW-2kV) |0,0,n) 
A2iip^'P^^'''^^)\l,l,n) = (2fi;2(r2-l) cos((^o) + 2fi:W + 2^;^) |1, l,n) 

A2(<^("P^'*"'))|mixed ±,n) = (57) 

2fi;^(r^cr + (r^ + 1) cos((/9o) ± y^r"^ — Ar'^ sin^((y9o)) [mixed ±,n) , 

where we have introduced the following notation: 



e = ^, (58) 

r^ — l 

di = cos{ipi) e {-1,1} . 

From these eigenvalues one can apply the previously described procedure 
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and obtain the formulae for the masses M(n_i, no, ni, . . . , n^-i): 

cosh(M(0, 0, n)) = |S + re| (59) 

cosh(M(l, l,n)) = \E- re\ (60) 

cosh(M(mixed ±, n)) = ^ S ^ ^{6^ - 1){E'^ - 1) + r^e^ (61) 
with S = $2(0 -r^ecr , 

where $2(0 is defined in Eq. (|52D and the remaining symbols in (^81). 

Now let us discuss these results. All the dependence on the occupation 
numbers rii lies in the quantity r. On the other hand, both S and r depend 
on the choice ifi = 0, tt. From the definition of r one sees that its maximum 
positive value is obtained whenever all states having -t^i = 1 are empty and 
those with i)i = —1 occupied. The maximum negative value is attained in the 
opposite situation. In both cases, the maximal absolute value is the same: 
d — 1. Let us now comment briefly on the main features of the spectrum 
formulae: 

n_i = no sectors 

This sector contains the scalar state corresponding to n_i = uq = Ui = 0, 
whose mass is given by: 

cosh(M5) = |<l>2(0-o^l • (62) 

For |r| > 1 (^ > 0) the scalar ground state corresponds to zero momentum 
0=0. For |r| < 1 (^ < 0) it corresponds to the doubler state ipi = ii. The 
state with minimum mass within these sectors corresponds to \t\ = d — 1 
and (pi = 0, which is the state associated to the matrix 70. Its corresponding 
mass is: 

cosh(My) = |<l>2(0 - ^(ci - 1)1 . (63) 
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n_i + no = 1 sector 

Let us consider first even space-time dimensions. In that case this sector has 
hghter states than the previous one. Notice that the argument of the square 
root in Eq. (|l]) can be written as {9E - 1)^ - ((S - 9^ - r^e^). Thus, the 
lightest state corresponds to the maximum value of |r| = rf— 1. Furthermore, 
one can prove that the mass decreases with |S|, and hence the hghtest state 
corresponds to (fi = 0. Combining this with the maximal r one sees that this 
lightest state corresponds to 7 and 707 ( 75 and 7075 in 4 dimensions). Its 
mass is given by formula ( |6lD with S = $2(0 — {d — l)r^e and |r| = d — 1. 
The corresponding critical line, where this lightest meson becomes massless, 
marks the edge of the physical region. It also marks the boundary of validity 
of our formulae. The equation for this critical line is given by: 

^2{0 = d9 . (64) 

Along this line, the other meson states have a positive definite mass. The only 
exception occurs for ^ = 1 (e = 0). The latter is an interesting region, where 
all the meson masses are degenerate and dependent only on C,- All states 
become massless at ^ = |. Notice that the region corresponds to the limit 
r ^ 00, K — i> with nr fixed. If we look at the expression of the action in 
this limit we see that the normal Dirac term becomes negligible with respect 
to the Wilson term. We have then essentially the masses corresponding to 
the scalar-gauge theory. Unphysical features could be expected to arise in 
this limit, because, as we commented in the previous section, |r| > 1 would 
be forbidden by the requirement of Osterwalder-Schrader positivity at the 
lattice level. However, it can be checked that the meson propagators in this 
limit satisfy the positivity condition. 

If we solve in this even- dimensional case for the critical hopping k^ which 
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defines the critical line as a function of r we get, in particular, the result 
/T;^(r = 1) = l/4(i, in agreement with the well-known fi;^(r = 1) = 1/16 for 



For odd space dimensions the previous analysis has to be modified. Now, 
it is not possible to satisfy |r| = d— 1 and ipi = simultaneously, due to the 
requirement that the number of gamma-fermions must be even. Thus, the 
lightest state has a higher mass. This state is a mixed vector-tensor state (7, 
and 7o7j). The critical line is given by |$2(0 — {d — 1)9\ = 1 and corresponds 
to a massless vector meson. For r > 1 this massless state is obtained for 
0=0 and for r < 1 it corresponds to one direction being cpi = n. There are 
no other massless particles except for r ^ 00. 

After this general presentation, we will now restrict ourselves to the most 
interesting cases of 3 and 4 dimensions, where the continuum theory is su- 
persymmetric. We will start by considering d = 4. The formulae for the 
propagators and masses follow from the general formulae obtained previ- 
ously. The normalisation of our operators in this case is given in Table 1. As 
mentioned previously, if we consider correlation functions of ultralocal oper- 
ators the only non-zero ones for Majorana fermions are the scalar (S) and 
the pseudoscalar-axial (PA) sectors. For ease of access to the results we now 
collect the expressions for the scalar mass Ms and the lightest pseudoscalar 
mass Mp: 

cosh{Ms) = |$2(0±3| (65) 

cosh (Mp) = 9E- v/(^2_i)(22_i) + 9e2 ^ 

where the ± in cosh(M5') is for |r| < 1 and \r\ > 1, respectively. We recall 
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the expression for $2 and S in this dimension: 

Sx^ - 14x + 7 



2x(l-x) 

$2(0 -3A 



$2(X) 



^, e and 6^ are, as usual, the ones in Eqs. 

The critical line where the pseudoscalar particle becomes massless is given 
hj K = K,c{i^), where the critical hopping parameter k^ is given by: 



23r2 + 9 + 3V9r4 + 46r2 + 9 
896r^ 



«^c(0 = ^7^^r:i • (66) 



Our formulae are equal for fields in the fundamental and the adjoint 
representation of the group, so that up to a normalisation factor, they can be 



directly compared with the results obtained previously on the literature [16 
||19|| . Once the appropriate normalisation of the fields is taken into account, 
we find perfect agreement with the results obtained by previous authors. 

Now we discuss the d = 3 case, for which no explicit results have been 
given previously. The expression for the propagator does not follow from 
our previous formulae. In this case the basis of meson operators is given by 
the set of 2 by 2 matrices Si = {1, cij}, where (Xi are the Pauli matrices. 
To comply with the conventional indexing of the Pauli matrices, we will 
adopt the names (yji, (y92, V^s) for the lattice momentum coordinates, time 
being assigned by convention to coordinate 1. It is easy to see that here 
the propagator splits into two separate blocks: the scalar and the vectorial 
sector, corresponding to the unit matrix and Pauli matrices respectively. By 
inverting the matrix 02(^1 — A2{<f) we find the following results for the 
scalar Gssi'^) and vectorial Gyvif) propagators: 

Gssi^) = ^.,, ^? T^ (67) 
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Gvvif) - , ^3 52- 



where the matrix elements of M on the ctj basis are given by 



1 ^ 
j,k=l 

3 

M^- = M-J^e^ifcafc/^fc, (68) 

fc=i 

eijk is the completely antisymmetric tensor and the functions cij jSj are given 
by 

3 

1=1 

Pi = 2tr sm{ipi) . 

For c? = 3, the functions iif and $2 adopt the form: 

'5\ 4x - 5 



H{x) = -T] D 



R 



2 J x(a;-l)(6a;-5) 



*^(^' = 2.(1-.) • (™' 

Let us now consider meson masses. The general formulae (|59| ) -(|6l D are 
valid for this case. The restriction to gamma-fermion states which are eigen- 
states of the charge conjugation operator C simply eliminates the double 
degeneracy of all levels. Taking into account the necessary evenness of the 
number of gamma- fermions, we arrive at: 

cosh(Ms) = |<l>2-a| =1 S + crel (71) 

cosh(Mii) = \^2- cr9\ =\E-ae\ (72) 

cosh(M±J = ^H±v/(S2-l)(02_i) + ,2(4_^2) ^ (73) 
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where, as usual, a = J2i=2'^i^ ^^^ '^i — cos(v9j) takes values {—1,1}. H is 
defined below Eq. (^) . The analysis of the lightest meson follows the general 
presentation done previously for odd space-time dimension. 

Finally, we will consider correlation functions of the 3-gluino operators 
given in Eq. (PB]). We recall that depending on the choice of the invariant 
colour tensor (/ or d) the spin matrix Si can be chosen completely symmetric 
or completely antisymmetric under the exchange of the 3 spinorial indices. 
The general expression for the propagator was given in Section 2 (Eq. (|35|) ), 
and involves the inversion of the matrix 63(^)1 — A3. Although performing 
this inversion explicitly in general seems a fairly complicated problem, we 
will be able to perform it inversion for the simplest case. This corresponds 
to a completely antisymmetric Si matrix. Actually, for rf = 4 it gives the 
full antisymmetric subspace. We proceed by introducing the antisymmetric 
matrices A a, where A is an spinorial index which can be looked at as the 
spin components of a composite fermion (a spin 1/2 fermion in 4 dimensions). 
The form of A a is given by: 

('^ \AiA2^3 _ (jAiA2(jA3A _ (jAiA-j,qA2A , (jA2MqAiA f^^\ 

The advantage of these matrices is that the states Cg^l^^) are the basis of 
an invariant subspace under the action of A3. This allows us to perform the 
appropriate inversion in this subspace. If we label by G a^i the momentum 
space propagator between the states Si = Aa and Sj = Aa', we can write: 

G'fl' = -i?3(0(6-3Tr(I,p,„))x 

^[(63(0 - 2K'r{r' - l)a)I + 2tK'ir^ - 1) J^ sin(^^)7,r M (75) 

where Tr{Ispin) is the dimension of the spin space (equal to 2'^'^ for even d) 
and a = ^^^coslipfj^). The inversion of the matrix contained in the previous 
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formula can be done in the standard way for a fermion propagator: (a + 

Setting as usual the space-like momenta to or tt, we might extract the 
mass of this fermion state from the position of the pole in the propagator. 
The result is: 



cosh(Mi)= rSg ± V(H3)2 - e , (76) 

where we recall that e = l/(r^ — 1) and the symbol S3 stands for: 



--3 



- rea (77) 



with a = ^j cos{(f\^^^""' ). It can be shown that, although the mass formula 
depends on the sign of k and r, the spectrum does not. Without loss of 
generality, we can set k > 0; then, the lowest mass is obtained ioi (f = 
when r>lor— l<r<0, and for (fi = TT\fi when r < — 1 or < r < 1. 

Obtaining the remaining masses of 3 gluino states analytically in general 
is a complicated problem. Nevertheless, we do not need them for the phys- 
ical conclusions of this paper. Indeed as we saw in the case of mesons the 
only point in parameter space were the scalar and the pseudoscalar became 
degenerate in mass occurred for k ^ 0, r ^ 00 with nr = vA fixed. Ex- 
tracting the masses and propagators for all p-gluino states in this limit is 
fairly simple. As in the case of mesons all spin states are degenerate, and 
the masses only depend on a. By looking at the expressions (^,^) for the 
matrix to be inverted when considering a p-gluino operator, we immediately 
see that: 

Gpl - Ap "^^ {Qp - 2KPrPa)I . (78) 

Now, the eigenvalues, whose zeroes gives the masses, are explicit. Their value 
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is given by: 

co.h(M.)^.(( '^''-f-^) )'")-.. (79) 

where we have used the definition (|25|) of C, and introduced the function 

,(.)4(.+ (!izi)) . (SO) 

As p increases ( - — ~ ^ ~ ' j , the argument of the function $ in (^), in- 
creases. It is easy to see that within the relevant range (x > 1) the function 
$(x) is monotonously increasing. This allows us to prove the following in- 
equality: 

Mp > Mq for p> q (81) 

As a consequence for any value of ^ in the physical range [0, |] and ipi = 0, 
the mass of the 3-gluino states is higher than the mass of the 2-gluino states. 
In particular at the critical point ^ = |, the 3-gluino states are massive. 

It is also interesting to know whether there are points within the physical 
region (bounded by the 2-gluino critical lines) where some 3-gluino states be- 
come massless. For example, in Ref. [0 it is predicted that if chiral symmetry 
is broken spontaneously then so is supersymmetry, and a massless goldstino 
particle appears. We have explicitly checked that this does not occur within 
our framework for the 3 and 4-dimensional cases. The 3-gluino states are 
always massive. 

4 Conclusions and future outlook 

In this section we will summarise our results and discuss their implications. 
We have obtained the spectrum of N=l SUSY Yang-Mills on the lattice at 
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large number of colours N^ and strong coupling, by considering the hopping 
parameter expansion as a sum over lattice paths (random walks). We have 
resummed the expressions in the hopping parameter in a certain region en- 
closing the origin, for an arbitrary value of the Wilson parameter r. We 
have worked at zeroth order in (3, the pure gauge coupling constant. How- 
ever, Wilson's action for the gauge part can be added either as a trace in 
the fundamental representation or as a trace in the adjoint one, with the 
corresponding couplings related to match the same naive continuum limit. 
Indeed, if we choose the adjoint version of Wilson action, our results (prop- 
agators and masses) are valid to all orders in (3 (probably only within some 
region encircling the origin). This can be proven in the same way that one 
sees that the quenched approximation is exact in our case. 

We have given formulae for the propagators and masses of 2 and 3 gluino 
states. The 2-gluino masses do coincide with the results for the meson spec- 
trum in ordinary lattice QCD at strong coupling [^-[^, and obtained by 



means of the effective potential method. Our method is based on the ran- 



dom walk resummation technique |2^, |23|. This generalises the method of 



Kawamoto [|TB| for r 7^ 1. Both methods have their relative advantages and 
disadvantages and occasionally there have been some conflicting conclusions 
(See the U(l) problem discussion in Refs. |1^ and [^). The random walk 



method does not rely on specific assumptions about the symmetries of the 
saddle point solution. It is rather based on the resummation of convergent 
series. Convergence is simple to see, since the number of lattice paths of 
length L grows like a power of L and the matrices involved are bounded. 
The radius of convergence is given by the distance to the nearest singularity. 
In our case, we have two resummations involved. One on the pure spike 
paths, which converges provided k^ < M2d-r)(r^-i) ^ ^^"^ ^ second one whose 
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border line in four dimensions, is given by the critical line (Eq. (^)) where 
the pseudoscalar becomes massless. Furthermore, in our case, the replace- 
ment of the Pfaffian by the square root of the determinant can be rigorously 
justified if |k| < 2d(\r\+i') ' '^'^^^ ^^^^^ ^ ^^^ ^°'-' short compared to the other lim- 
its. Finally, we want to stress that we have provided formulae for arbitrary 
space-time dimensions d, which could be of interest for other researchers in 
the field. 

Apart from the technical interest of our methods and results, we consider 
that the main usefulness of our results, is that they provide a guideline and 
boundary formulae for groups investigating this model by Monte Carlo sim- 
ulations. Of course, our results are only valid for large N^.-, but experience 
teaches us that this is frequently a numerically good approximation. There 
is one issue in which unfortunately our method could perhaps not help. It 
has been predicted, that this model should have a first order phase transition 
in K |]2^ . Our work predicts the presence of a second order phase transition 
associated to the vanishing of the pseudoscalar mass. It can be argued how- 
ever, that a series expansion like ours can be seen as the expansion around 
one of the effective potential vacua. Thus, the mentioned second order tran- 
sition could lie in the metastable phase. The point about the order of the 
phase transition should be settled by future Monte Carlo simulations. 

Finally, it is tempting to speculate about the relevance of our results in 
the spirit of supersymmetry breaking. For that purpose one is interested 
in critical lines where a continuum limit can be defined. The states whose 
lattice mass vanish at the critical line, are the states that survive this con- 
tinuum limit. If supersymmetry is recovered at this continuum limit one 
expects these states to form a supermultiplet. The analysis of Curci and 
Veneziano leads to a multiplet in which in addition to the pseudoscalar par- 
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tide, there is a scalar one and a spin | fermion. These particles have equal 
positive continuum masses. Since the contribution to the mass of the pseu- 
doscalar comes through the anomaly, which vanishes in our case, we should 
expect a massless multiplet. Nevertheless, it is doubtful that the analysis 
of Curci and Veneziano applies at strong coupling since it is based on the 
continuum SUSY Yang-Mills lagrangian. By power counting, this model has 
a single relevant parameter, the gluino mass, and hence fine tuning one of 
the bare couplings one could find a line corresponding to vanishing mass and 
restored supersymmetry. However, at strong coupling the gluonic degrees of 
freedom stay of the order of the cut-off. Hence, the low energy lagrangian, if 
super symmetric, would rather coincide with the Wess-Zumino model. This 
has many more relevant operators (the different masses and couplings) and 
demands tuning of more bare parameters to approach it. In this respect 
the situation in 3 dimensions could be interesting since the model would be 
interacting. In 4 dimensions we would expect a free low energy lagrangian 
giving the physics of the continuum limit. With all these concerns in mind 
we did not want to loose the opportunity to explore the k — r parameter 
space in search for degenerate low energy multiplets. Actually, we concluded 
that the only point where several mesons become massless is in the limit 
K — > 0, r — * oo and rn = J-prr (i-^-, ^ = |). The masses at this point are 
the ones corresponding to a gauge-Higgs system: an interesting model in its 
own right. At this point we looked at the masses of the p-gluino states with 
p > 2. This includes fermionic degrees of freedom (for p odd). However, 
we showed that these states remain massive (in cut-off units) at this critical 
point. So that the main conclusion is that there is no point in the /t — r plane 
giving a possible candidate for a supersymmetric continuum limit. 

We conclude the paper by mentioning a few possible improvements of our 
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results. First of all, the possibility of extending the calculations and formulae 
to next to leading order in l/N^ seems a feasible one. The most important 
consequence of this extension could be in cases when the effects are absent 
to leading order, like the effect of quenching, anomalies, etc. Then one can 
try to include higher orders in (3 fundamental, or combined l/N^ and Padjoint- 
Then, of course, it would be very good to rederive the results of this paper 
with the effective action method. This technique, as mentioned previously, 
would allow the discovery of possible first order phase transitions. Finally, it 
should be commented that our methods and results could be used to study 
other supersymmetric models, such as SUSY QCD. 
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Appendix 

In this appendix we will present the terminology and main results on random 
walks that we will need in the text. Not to conflict with other definitions, 
we will refer to these random walks as lattice paths and a precise definition 
will be provided. We will work in arbitrary dimension d and will employ 
additional definitions given at the beginning of Section 2. Proofs will not be 



given. For that we address the reader to Ref. |]2^ . 
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A lattice path of length L is an element 7 = (n, a) G C ^ I^ . The 
point n G £ is the origin of the path, a is the path sequence, and m = 
n + X]i=i "l^(Q^i) is the endpoint of the path. We can label the different path 
sets as follows. Let Siin) be the space of all paths with origin n and length 
L] SlIu — > m) is the space of all paths with origin n, endpoint m and length 
L. Now we will introduce the notion of a spike. 

A path 7 = (n, «!,... , ai) has spikes if there exist one integer 1 < i < L 
such that ctj+i = di. In the converse case one says that the path has no 
spikes. The set of all paths without spikes of length L and origin n is 
labelled Siin) [Siin —>■ m) if the endpoint is fixed to m). Now, we can 
associate to any path 7 a corresponding path 7 called its reduced path, by 
simply eliminating in an orderly manner all pairs Oj+i = di in its sequence. 
Notice that this procedure preserves the origin and endpoint of the path. 

Our main strategy will be to convert the sum of paths into a sum over 
reduced paths. For that we would first need to compute the following func- 
tion: 

00 

F{L,z)=J2^'N{L,p) , (82) 

where N{L,p) is the number of paths of length L + 2p whose reduced path 
is a path (with no spikes) of length L. This number does not depend on the 
path itself but only on its length. In Ref. P3[ it is shown that: 



^'"-'^ (i-kmi-V '^'' 



^ ^ i-ViziPl^ilf . (84) 

In addition, we would need to be able to perform sums over the set of 
reduced paths with origin in n of a product of matrices. To be specific, let 
Aa for a G /, be a collection of matrices satisfying A^Ao = AI, with I the 
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unit matrix. Then we want to compute the matrix: 



r(A) = l + 5^ Y, A,, ■••A,,. (85) 

L=l {n,a)e5L(n) 



It can be shown [^] that T(A) is given by the formula: 

r(A) = (l-A)(l + (2rf-l)A-A)-i , (86) 

where A = Y^^^j A^. 

In the text we will need to evaluate an expression like Eq. (|85[) with a 
slight variation. We would need to sum only over paths going from one lattice 
point X to another one y. This modified situation can be reduced to the case 
given before by the following procedure. Instead of considering the matrix 
Aq, we will multiply it by a phase e**^", where cpp, = —f^. Then, diagrams 
that go from x to y have coefficients that go like e*'^'-^"^-' . Thus, the required 
expression can be obtained by projecting onto this term: 

%^y{A) = ll(^J ^^ e-^^^-yHl-X){l + {2d-l)X-A'{<p)r' , (87) 
with A'(0) = J2a^j e'^- A„. 
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Table 1: Nomenclature and normalisation for the spin Si and colour Cj matri- 
ces used to define the four dimensional 2-gluino operators. The last column 
gives the number of components associated to each one of them. The vec- 
tor and tensor spin matrices V{p) and T{pa) are symmetric, thus forbidden 
when the field is Majorana. 
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